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Abstract 

Starting from the developed generalized point process model of 1// noise (B. Kaulakys et al, 
Phys. Rev. E 71 (2005) 051105; cond- mat/0504025 ) we derive the nonlinear stochastic differential 



equations for the signal exhibiting 1/ f@ noise and 1 /x A distribution density of the signal intensity 
with different values of (3 and A. The processes with l/f@ are demonstrated by the numerical 
solution of the derived equations with the appropriate restriction of the diffusion of the signal 
in some finite interval. The proposed consideration may be used for modeling and analysis of 
stochastic processes in different systems with the power-law distributions, long-range memory or 
with the elements of self-organization. 
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I. INTRODUCTION 



The entirely uncorrelated in time stochastic signals exhibit white noise - the power spec- 
tral density S(f) as a function of the frequency / is constant, while the Brownian motion 
of the signal intensity without correlations between increments results in 1/f 2 and in the 
Lorentzian power spectra. The widespread occurring signals and processes with 1/f spec- 
trum (see, e.g. □DDDD and references therein) cannot be understood and modeled in 
such a way. 

"1// noise" is a type of stochastic processes which power spectral density at low fre- 
quencies behaves like S(f) ~ 1//^, where the exponent f3 is close to 1. In contrast to the 
Brownian motion and 1/f 2 noise generated by the linear stochastic equation, simple systems 
of differential equations, even linear stochastic equations, generating signals with 1/f noise 
are not known. 

Recently, starting from the simple point process model of 1/f noise flEl, we derived the 
nonlinear stochastic differential equation 

^ = x 4 + x%) (1) 

for the signal intensity x generating processes with pure 1/f noise and the inverse cubic, 
P(x) ~ 1/x 3 , distribution of the signal intensity Here t s is the scaled time and £(t s ) is 
a Gaussian white noise satisfying the standard condition 

(mm) =*(*.-o (2) 

with the brackets (...) denoting the averaging over the realizations of the process. 

The aims of this paper are the derivation of the class of the stochastic nonlinear differential 
equations exhibiting 1/f 13 noise and l/x x distribution density of the signal with different 
values of the exponents (3 and A and the numerical demonstration of the proposed model 
for generation of the long-range fractal processes. 

II. THE MODEL 

We start from the point process as a sequence of correlated pulses or series of events 

I(t) = aJ2^(t-t k ). (3) 



Here S(t) is the Dirac ^-function and a is a contribution to the signal or current of one pulse 
at the time moment tk. We consider the stochastic multiplicative process for the interevent 
time T k = t k+1 - t k (a |<| 

Tfc+i = T k + JT^' 1 + O-T-fe^fc, (4) 

where the (average) interevent time fluctuates due to the random perturbation by a sequence 
of uncorrelated normally distributed random variables {e k } with zero expectation and unit 
variance, o denotes the standard deviation of this white noise and 7 <C 1 is a coefficient of the 
nonlinear damping. It has been shown analytically and numerically that the process 

(0) and (jlj may generate signals with the power-law distributions of the signal intensity and 
1/f 13 noise. 

Transformation of Eq. (jlj) to the Ito stochastic differential equation in fc-space is 

Transition from the occurrence number k to the actual time t according to the relation 
dt = r k dk yields 

^ = 1 r^ + ar^(t). (6) 
The standard transformation 10] of the variable from r k to the averaged over the time 

n 

interval r k intensity of the signal x = a/r k by analogy with Ref. [£| yields the stochastic 

nonlinear differential equation 

fj r r 4-2^t 5/2-/1 



Introducing the scaled time 
and the new parameters 



a 2 

,3-2^' 



(8) 



1 = *-* F = 1 -^ ( 9 ) 



5 r 1 7 

--11, r = 1 - - 

I a 
we obtain the class of Ito stochastic differential equations 

^ = Yx 2 ^ 1 + x^{t s ). (10) 

Eqs. (jlUJl . as far as it corresponds to the analyzed in Refs. point process ©-(01), should 
generate the signals with the power-law distributions of the signal intensity, 

P(*)~-i, A = 2fo-r), (11) 

Jb 
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FIG. 1: Typical examples of the solutions of Eqs. (|14jl and (j!5j) : a) with the parameters T = 0.75, 
rj = 2, and k = 0.1 and b) with the parameters T = —0.2, r] = 1.5, and k = 0.01. 



and 1/ f 13 noise, 



1 2F 4- 1 



(12) 



" 277 -2 

According to the general theory [lfj the exponentially restricted diffusion with the dis- 
tribution densities 



P(x) 



x/ 



exp 



X 



X 



should generate the stochastic differential equations 



dx m ( x™: 



x 



m—l+2ri 



dt, 2 \x m+1 ~ 2r i xl 



+ fx 2 ^ 1 + x^{t s 



where m is some parameter. 



III. NUMERICAL ANALYSIS 



(13) 



(14) 



For the numerical solution of Eq. (|14|) we can take the integration steps from the equation 
x^y/Ki = KXi, with k < 1 being a small parameter. This corresponds to the case when the 
change of the variable Xi in one step is proportional to the value of the variable. As a result, 
we obtain the system of equations 



2 I XT 



X; 



x : 



t+1 



U + 



X 



K 

277—2 ' 



(15) 



4 




FIG. 2: a) Numerically simulated distribution density of the variable x according to Eqs. (|15j) . open 
circles, compared with the expected distribution density solid line; and b) power spectral 

density, obtained from the numerical solution of Eq. (|15|) . open circles. Solid line in b) represents 
the analytical power spectral density slope (|T2*|). Parameters used are x m i n = 1, x max = 10 3 , m = 1, 
r = 0.75, r\ = 2, and k = 0.1. The signal was calculated from N x = 10 6 points and averaged over 
n = 100 realizations, c) and d) represent the distribution density and the power spectral density, 
respectively, with the parameters m = 1, T = —0.2, r\ = 1.5, and k = 0.01. 

In ngureQthe typical examples of the signals as solutions of Eqs. (114)) and ()15|) are shown. 
The distribution densities P(x) of the variable x, obtained by the numerical simulation of 
Eq. ()15|). are shown in figures|21a) and c). The power spectral densities S(f) are shown in fig- 
ures|21b) and d). Numerical simulation of distribution densities and power spectral densities 
are in good agreement with approximate expressions, Eq. (fTTjl and Eq. (fT2j) . respectively. 
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IV. CONCLUSIONS 



We derived and analyzed a class of stochastic nonlinear differential equations for the 
signal exhibiting 1/ f 13 noise and l/x x distribution density of the signal in any desirable 
wide range of frequency and of the signal intensity. The proposed technique may be used for 



modeling of the stochastic processes in different systems (e.g., in financial systems |9 



and the Internet 



11 



with the power-law statistics of the signal characteristics. 
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